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Unlike our basic theories of space and time, quantum mechanics is not a locally causal theory. 
Moreover, it is widely believed that any hopes of restoring local causality within a realistic theory 
have been undermined by Bell's theorem and its supporting experiments. By contrast, we provide 
a strictly local, deterministic, and realistic explanation for the correlations observed in two such 
supporting experiments, performed independently at Orsay and Innsbruck. To this end, a pair of 
local variables is constructed to simulate detections of photon polarizations at various angles, chosen 
freely by Alice and Bob. These generate purely random outcomes, srf = ± 1 and SB = ± 1, occurring 
within a parallelized 3-sphere. These outcomes do not depend on the contexts of measurements, but 
are determined entirely by the randomly chosen initial orientation of the 3-sphere. When they are 
compared, however, the correlation between them turn out to be exactly equal to — cos2(a — j3), 
with the corresponding CHSH inequality violated for the polarization angles a, a', f3, and /3' in 
precisely the manner predicted by quantum mechanics. The key ingredient in our explanation is the 
topology of the 3-sphere, which remains closed under multiplication, thus preserving the locality 
condition of Bell. It allows us to model the physical space as a 3-sphere, and reveals that the illusion 
of quantum nonlocality in the present case stems from a twist in the Hopf fibration of the 3-sphere. 

(N _ 

i— | ' One of the first steps we often take towards measuring a physical quantity is to set up a Cartesian coordinate system 
i-C . {x,y,z} in the Euclidean space E3. This amounts to modeling the Euclidean space as a 3-fold product of the real 
line, IR 3 . This procedure has become so familiar to us that in practice we often identify E3 with its Cartesian model, 
and simply think of IR 3 as the Euclidean space. As we shall see, however, this seemingly innocuous act of convenience 
comes with a heavy price: It is largely responsible for the illusion of quantum nonlocality [HQ- I n what follows we 
^ ■ shall demonstrate that, once a correct coordinate- free geometric model of the physical space is used, the correlations 
, observed in the EPR-inspired experiments involving photon pairs Q ^t] Q , namely 
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sf(a) = ±1, 3&{fi) = ±1, 



> ■ 8(a) = 0, £{p) = 0, 

^f- 8{a,p) = -cos2(a-/3), (1) 



are easily understood, in a strictly local-realistic terms. Here £ (a, 0) represents the expected value of simultaneously 
observing remote measurement results si (a) and 3$(/3) along the polarization angles a and (3, respectively. 



Euclid himself of course did not think of E3 in terms of triples of real numbers. He defined its representation 
axiomatically, entirely in terms of primitive geometric objects such as points and lines, together with a list of their 
. . ■ properties, from which his theorems of geometry follow. Today we know, however, that it is quite problematic to 
give a suitable definition of Euclidean space in the spirit of Euclid, and hence in physics we instinctively identify E3 
with IR 3 whenever possible. But there is no natural, geometrically-determined to identify the two spaces without 
;_h , introducing an unphysical notion of arbitrarily distinguished coordinate system [6j. This difficulty is clearly relevant 
in the study of Bell's theorem [2[, for time and again we have learned that surreptitious introduction of unphysical 
ideas in physics could lead to distorted views of the physical reality. A coordinate-free representation of the Euclidean 
space is undoubtedly preferable, if what is at stake is the very nature of the physical reality. 

Fortunately, precisely such a representation of E3 , with a rich algebraic structure, was provided by Grassmann in 
1844 0. As in Euclid's geometry, the basic elements of this powerful structure are not coordinate systems, but points, 
lines, planes, and volumes, all treated on equal footing. Today one begins this framework by postulating a unit volume 
element (or a trivector) in E3, defined by 

I = e x A e y A e z , (2) 

with {e^, e y , e z } being a set of orthonormal vectors, "A" the outer product of Grassmann, and I 2 = —1 [7]. Each 
vector ej is then a solution of the equation I A ej = 0, and every pair of them respects the fundamental geometric 
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product 

ej e fc = ej ■ e k + e 3 A e k , (3) 

where " • " represents the inner product, and ej A are unit bivectors (which also square to —1) with counterclockwise 
sense for the cyclicly permuted indices (j, k — x, y, or z). The resulting real geometric structure is a linear vector 
space spanned by the basis 

{1, G x , e.y, G z: G x A Gy : Gy A 6 2 , G 2 A Q X 7 f\ A ^z}; (4) 

which encodes a graded linear algebra of dimensions eight. This algebra intrinsically characterizes the space E3. 

Our interest, however, lies in a certain subalgebra of this algebra, the so-called even subalgebra of dimensions four, 
defined by the bivector (or spinor) basis 

{1, e x A e y , e y A e z , e z A e x }. (5) 

Crucially for our purposes, this subalgebra happens to remain closed under multiplication. Consequently, it can be 
used by itself to model the physical space. In fact, it provides the most natural coordinate-free representation of the 
physical space as a parallelized 3-sphere [6[ , which differs from IR 3 only by a sing le point: S 3 = IR 3 U {00}. Both H 3 
and S 3 are three-dimensional manifolds that are parallelized by a vanishing Riemann curvature, but unlike in IR 3 the 
torsion within S 3 is non- vanishing Q . S 3 is thus a one-point compactification of H 3 that remains as flat as H 3 Q . 

The vectors and trivectors are no longer intrinsic to the above subalgebra, but belong to a dual space. Only the 
scalars and bivectors — treated on equal footing — are taken to be intrinsic parts of the subalgebra. This can be seen 
more clearly if we use the condition / A e ^ = to rewrite the basis bivectors defined in Eq. ([5]) as I ■ e z , I ■ e x , and 
I ■ e y . Their geometric product, analogous to the one in Eq. ([3j) , then leads to the defining equation of this subalgebra: 

(/ ■ e,-) (/ ■ e fe ) = - 6 jk - e jk i {I ■ e t ). (6) 

Evidently, despite the occurrences of trivectors and basis vectors, only the basis scalar and bivectors are involved in 
this definition. Consequently, in what follows only scalars and bivectors (and their combinations) will have direct 
physical significance — vectors and trivectors will merely facilitate computational ease, or "hidden" variables. 

Given the bivector basis defined by Eq. ([5]) , any generic bivector such as / • a can be expanded in this basis as 

/ ■ a = { a x e y A e 2 + a y e z A e x + a z e x A e y }. (7) 

In many ways bivectors — and also vectors and trivectors — behave just like ordinary numbers within this geometrical 
framework. For this reason Grassmann referred to them as "extensive magnitudes" , and treated them with the same 
respect as he would treat ordinary numbers. And for the same reason Hestenes today refers to all such quantities as 
"directed numbers" 0]. They are simply numbers of higher grades and dimensions, with built-in directional attributes. 

It is also worth noting that, although there is isomorphism between the vector subspace and the bivector subspace, a 
bivector is an abstract entity of its own, with properties quite distinct from those of a vector Q ■ In fact it is a skeleton 
of only three properties: (1) a sense, quantified by a + or — sign, indicating whether it represents a counterclockwise 
or clockwise rotation, (2) a magnitude, which is equal to unity for all bivectors we will be considering, and (3) a 
direction, which can be specified by a dual vector normal to its shape-independent plane. Thus, despite appearances, 
neither the trivector / nor the vector a is an intrinsic part of the bivector / ■ a. Moreover, it is easy to verify that 
unit bivectors such as / • b simply represent intrinsic points of a unit 2-sphere contained within the 3-sphere: 

||/.b|| 2 = (+/.b)(-/.b) 

= -7 2 bb = bb = b- b = ||b|| 2 = + 1. (8) 

More precisely, they represent the equatorial points of a parallelized 3-sphere modeling the physical space Q . 
Given two such unit bivectors, say / • a and / • b, the bivector subalgebra ([6]) leads to the well known identity 



(/ • a)(I • b) = - a- b - I ■ (a x b), 



(9) 
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provided we use the duality relation a A b = J ■ (a x b) defined between the vector a x b and the bivector a A b. 

The above identity provides a natural representation of points of a parallelized 3-sphere [8| . The bivectors I ■ a and 
I ■ b appearing on its LHS represent the equatorial points of the 3-sphere, and the real quaternion appearing on its 
RHS — which is a sum of a scalar and a bivector — represents a non-equatorial point of the same sphere. An equator of 
a parallelized 3-sphere, however, which is a 2-sphere, does not remain closed under multiplication. The 3-sphere itself, 
on the other hand, does remain closed under multiplication, thereby correctly encoding the topology underlying the 
subalgebra ((6]) . Conversely, any arbitrary point P of a parallelized 3-sphere can always be factorized into any number of 
points: P = ABCD... Needless to say, this is a highly nontrivial and powerful property of the 3-sphere. To appreciate 
its non-triviality, consider a product of infinitely many points of a 3-sphere. Such a product will simply be another 
point of the 3-sphere By contrast, this will not be true in the case of a 2-sphere even for just two points, as is 
evident from the above identity. The fact that both 3-sphere and its algebraic representation ([9J remain closed under 
multiplication suggests that we should represent measurement results in the present case by local maps of the form 

j*(a, A) : H 3 x A — > S 2 C S 3 , (10) 

where A is a space of complete states A, vector a 6 H 3 specifies the context of measurement, and S 2 is an equatorial 
2-sphere within a unit, parallelized 3-sphere. Then not only the measurement results, but also their products would 
remain within the 3-sphere, thereby respecting the locality (or factorizability) condition of Bell [2] . We shall see that 
a local-realistic interplay between the points of such a 3-sphere and its equatorial 2-sphere is what is truly responsible 
for the EPR correlations manifested in nature. 

So far we have considered the bivector subalgebra (j6]) with arbitrarily fixed basis, as in definition ([5]). The convention 
usually is to assume a right-handed set of basis bivectors, and so far we have followed this convention. The algebra 
itself, however, does not fix the handedness of the basis. We could have equally well started out with a left-handed 
set of bivectors, by letting — / instead of -I- 1 fix the basis. Equation ([6]) would have then had the alternate form: 

(-1 ■ Bj) {-I ■ &k) = - Sjk - ejki {-I ■ ei). (11) 

Comparing this equation with equation (jB]) we see that there remains a sign ambiguity in the definition of our 
subalgebra (cf. Refs. and Q): 

(/ • e,) (/ • e fc ) = -5 jk ± e jH (J-e,). (12) 

Consequently, following the time-honored mathematical practice of turning an ambiguity of sign into virtue, we define 
the handedness of this entire subalgebra as our "hidden variable." In other words, we specify the complete state of 
the photons we are about to study as — ±1, thereby defining the basis of our entire subalgebra by the equation 

(/x • ej) (/x • efc) = - S jk - ejki (fi ■ ej). (13) 

The identity ((9|) for the generic bivectors then becomes 

( (J, ■ a)( /it • b) = a • b fj, ■ (a x b), (14) 

along with indefinite duality relation a A b = fi ■ (a x b). In other words, the duality between the wedge product 
and cross product within our subalgebra alternates between the right-hand and left-hand rules Q . Clearly, then, our 
complete state fi = ± I represents a far deeper hidden structure than the variables considered by Bell [2} ■ It un-fixes 
the orientation of the entire physical space S 13 , and turns it into a shared randomness between Alice and Bob. 

We are now well equipped to take up the question of EPR correlations exhibited by the pair of entangled photons in 
the actual experiments 0101 • This question, of course, has been well scrutinized in the literature We will restrict 
to the most basic aspects of the question, and follow its treatment given in Refs. Q and 0. In a quantum mechanical 
description of the experiment involving photon pairs one usually assumes that the system has been prepared in the 
singlet state 

|*_) = -^=-{\H) 1 ^\V) 2 ~ \V) 1 ^\H) 2 ], (15) 

where | H ) and | V ) denote the horizontal and vertical polarization states of the photons along the directions e x and 
e y , respectively; and the subscripts 1 and 2 refer to the photons 1 and 2, respectively. The photons are thus assumed 
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to be propagating in the e z direction. We could equally well consider the polarization state | but that would not 
add anything significant to our concerns here. Both polarization states, | and | 9—}, are invariant under rotations 
about the e z axis, but the state | is even under reflections, whereas the state | Vt-) is odd under reflections. 

In a typical experimental run Alice and Bob measure polarizations of the photons along two different directions 
in the plane perpendicular to the e z axis. Alice measures polarizations along the direction a, which makes an angle 
a with the e^, axis, whereas Bob measures polarizations along the direction b, which makes an angle /? with the e x 
axis. Individually, the binary results observed by Alice and Bob, namely A(a) — ± 1 and B(f3) = ± 1, are found to 
be completely random, with equal probabilities for the outcomes + 1 and — 1. When these results are compared, 
however, they are found to be strongly correlated, in agreement with the quantum mechanical predictions we have 
summarized in Eq. ([I]). 

Our goal now is to reproduce these quantum mechanical predictions exactly, within the geometrical model of the 
physical space discussed above. To this end, we have assumed that the complete state of the photons is given by 
fi = ±1, where / is the fundamental trivector defined in Eq. $2$. The detections of photon polarizations observed 
by Alice and Bob along their respective axes a and b, with the bivector basis fixed by the trivector /x, can then be 
represented intrinsically as points of the physical space S 3 , by the following two local variables: 

S 3 3^(a,») = (-/.a)( + M-5) = | + ! * fl = +I (16) 

1—1 if fx = —I 

and 

S 3 3^(/3,/x) = (+/•£)( + /*• b) = (~J t fl = + i 

|^+1 if (i = - I , 

with equal probabilities for /j being either + 1 or — I, and the rotating vectors a and b defined as 

a = cos 2a + e y sin 2a , 

and b = cos 2/3 + e y sin 2(3 . (18) 

Note that s?/(a, fi) and 0K/3, fi), in addition to being manifestly realistic, are strictly local variables 1 . In fact, they 
are not even contextual [10(. Alice's measurement result, although refers to a freely chosen angle a, depends only on 
the initial state fx; and likewise, Bob's measurement result, although refers to a freely chosen angle /?, depends only on 
the initial state fi . In other words, all possible measurement results at all possible angles are completely determined 
by the initial orientation of the physical space specified by /x, or equivalently by that of the 3-sphere. Moreover, it is 
easy to check using the identity ([Tip that srf (a, fi) and 3&[fi, \i) are purely binary numbers. For consider the product 

(+/• £)(+!• a') = - a- a' - I - (ax a') 

= — cos 6*55/ — (I -c) sin / , (19) 

where c = a x a'/|a x a'|, and 9^a' is the angle between a and a'. The RHS of this identity is a quaternion, which 
represents a non-equatorial point of the parallelized 3-sphere [8] . Physically, it represents a rotation by angle 29^ i 
about the c-axis. The limit a' — > a then gives (+ J • a) (+/ • a) = — 1, which shows that unit bivectors square to — 1. 
Moreover, in equation (j8]) we have already seen the conjugate limiting case; namely, (— I ■ a)(+ I ■ a) = + 1. These 
two cases show that the pairs of numbers { + 1, — 1 } are intrinsic parts of the parallelized 3-sphere. They are in fact 
antipodal points of one of its equators, representing the set of measurement results observed by our experimenters. 

In statistical terms, however, the above variables are raw scores as opposed to standard scores [ll[. Recall that a 
standard score indicates how many standard deviations the observation or datum is above or below the mean. If x is 
a raw (or unnormalized) score and x is its mean value, then the standard (or normalized) score, z(x), is defined by 

z (x) = ^, (20) 

cr(x) 



1 Needless to say, srf(a, fi) and 3%{fi, /x) are two different functions of the random variable /x. Moreover, they are statistically independent 
events occurring within a 3-sphere, with factorized joint probability P(s$ and £$) = P(£/)xP{£&) < ^. Therefore their product si 38 is 
guaranteed to be equal to —1 only for the special case a = (5. For all other a and 0, si 38 will alternate between the values —1 and +1. 
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where <r(x) is the standard deviation of x. A standard score thus represents the distance between a raw score and 
the population mean in the units of standard deviation, and allows one to make comparisons of raw scores that come 
from very different sources. In other words, the mean value of the standard score itself is always zero, with standard 
deviation unity. In terms of these concepts the bivariate correlation between raw scores x and y is defined as 



£(x, y) 



lim 

n»l 



1 " 



x)(y l -y) 



i = 1 



cr(x) cr(y) 



lim 



1 " 

-£z(xW) 

71 J 



(21) 
(22) 



It is vital to appreciate that covariance by itself — i.e., the numerator of the equation (|2ip by itself — does not provide 
the correct measure of association between the raw scores, not the least because it depends on different units and 
scales (or different scales of dispersion) that may have been used (advertently or inadvertently) in the measurements 
of such scores [llj |. Therefore, to arrive at the correct measure of association between the raw scores one must either 
use equation (|2ip . with the product of standard deviations in the denominator, or use covariance of the standardized 
variables, as in Eq. (|221) . 



These basic statistical concepts are crucial for understanding our local-realistic explanation of the EPR correlations. 
As we saw above, the variables .s/(a, fi) and ■$§?(/?, /x) are pure binary numbers, albeit occurring as antipodal points 
within a parallelized 3-sphere. As random variables, however, they are products of two factors — one random and 
another non-random. For instance, within srf (a, /z) the factor ( + fi ■ a ) is a random bivector — a function of the hidden 
variable /i, whereas the factor (—7-5) is a non-random bivector, independent of the hidden variable fi. Consequently, 
as a function, each of the numbers .s/(a, fi) and ^(/3, (J.) is generated with a different standard deviation — i.e., a 
different size of a typical error. More specifically, the number £/(a, fj,) is generated with the standard deviation 
(—/■a), whereas the number i^(/3, fi) is generated with the standard deviation (+ 1 ■ b). These deviations can be 
calculated easily. Since errors in a linear relation such as (fT6|) propagate linearly, the standard deviation of &/ (a, fi) 
is equal to (— I ■ a) times the standard deviation of ( + fi ■ a) (which we write as a(A)), and similarly the standard 
deviation of 38{J3, n) is equal to (+ 1 ■ b ) times the standard deviation of ( + /x • b ) (which we write as a(B)): 



and a(m) = ( + I-b)a(B). 



(23) 



But since all bivectors we have been considering are normalized to unity, and since the mean of ( + fi ■ a ) vanishes on 
the account of fi being a fair coin, its standard deviation is easy to calculate, and it turns out to be equal to unity: 



a(A) 



\ 



- V \\A(a, n*) - A(a, 



n 

x -Eik+^-)-°ii 2 =1. 



(24) 



where the last equality follows from the fact that all bivectors ( + fi 1 • a ) are normalized to unity. Similarly, we find 
that the standard deviation cr(B) is also equal to 1. As a result, the standard deviation of &/(a, fi) works out to be 
equal to (— J • a), and the standard deviation of ^(/3, £t) works out to be equal to (+ I ■ b). Putting these two results 
together, we arrive at the following standard scores corresponding to the raw scores (TIB)) and (jTTJ) : 



A(a, fi) 



(a, fj,) 



(a, n) 



£?(a, n) - 



(-/■a) 



( + /x-a) 



(25) 
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and 



B(/3, fx) 



m) - m> /*) 



ct(^) 

(+/.b) 



( + /x-b), 



(26) 



where we have used the identities (+ 1 ■ a)(— 7 • a) = +1 and (— 7 ■ b)(+ 7 ■ b) = +1, respectively (see Appendix for 
the formal basis of these results). 

Not surprisingly, just like the raw scores (a, fx) and 3$(0, fi), these standard scores are also strictly local variables: 
A{a 1 fx) depends only on the freely chosen angle a and the initial state fx, and likewise B(/3, fx) depends only on the 
freely chosen angle (3 and the initial state fx . Moreover, despite appearances, A(a, fx) and 77(/3, fx) are simply binary 
measurement results, ± 1, albeit occurring within the compact topology of the 3-sphere rather than the real line: 



S 3 D S 2 3 A(a, fx) = ( + ■ a) = ±1 about a e R 3 , 
S 3 D S 2 3 B(/3, fx) = ( + fx-b) = ±1 about b e K 3 . 



(27) 
(28) 



In fact, since the space of all bivectors ( + fx ■ a) is isomorphic to the equatorial 2-sphere contained within the 3-sphere 
each standard score A(a, fx) of Alice — corresponding to a predetermined element of reality of the photon — is 
uniquely identified with a definite point of this 2-sphere. In other words, every predetermined polarization of a given 
photon gives rise to a unique standard score A(a, fx), which in turn is unambiguously represented by a definite point 
of the equatorial 2-sphere, and likewise for the standard scores B(P, fx) of Bob 0. 

Now, since we have assumed that initially there was 50/50 chance between the right-handed and left-handed 
orientations of the physical space — i.e., equal chance between the initial states fx = + 7 and \i = — I, the expectation 
values of the local outcomes trivially work out to be 



£{6) 



lira 

n>l 



1 



H) = 0, 



(29) 



where 8 — a or j3. On the other hand, to determine the correct correlation between the joint observations of Alice 
and Bob we must calculate covariance between the standard scores A(a, fx) and B(/3, fx), not the raw scores: 



£ (a, j3) — lim 



E 



A(a, li 1 ) B{p, li 1 ) 



(30) 



Conversely, the covariance of the raw scores will not give us the correct correlation between the observations of Alice 
and Bob, because (a, li) and li) have been generated with different scales of dispersion for each direction, 

and hence they are subject to different standard deviations or random errors for each direction 11} . On the other 
hand, the standard scores — by definition — have been standardized for all kinds of possible random errors. Moreover, 
their product A(a, li) S(/3, li) — as a non-equatorial point of the 3-sphere — represents a definite binary value: 



S 3 3 A(a, li) B(j3, fx) = ± 1 about a direction in IR 4 



(31) 



In fact, our 3-sphere is entirely made of such product points, each of definite value + 1 or — 1, depending on its 
orientation defined by fx. More precisely, the space of all such product points P = AB is nomotopic to a round 
3-sphere: P x 2 + P| + P 3 2 



Pf = 1 



12j . This can be seen more clearly if we expand AB using the identity (fH)) : 



A(a, fx)B{(3, fx) = ( + fi-a)( + fx-b) 

= — a ■ b — fx ■ (a x b) 

= - cos2(a - p) + (fx ■ e z ) sin2(a - 0). 



(32) 



Evidently, the product AB describes a circle of points within S 3 [l^, each of definite value + 1 or — 1, depending on 
a, f3, and fx ■ e z e S 2 C S 3 . Substituting this parameterization of AB into Eq. (j3"0")l . the correlation between the raw 
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scores (a, fi) and ^(/3, /x) can now be easily seen to reproduce the corresponding quantum mechanical prediction: 



£(a, P) = - cos2(a - f$) + Km 

n » 1 

= -cos2(a-/3) + 0. 



n 

- V(^-e,) sin2(a-/3) 



i = 1 



(33) 



Here the summation over microstates fi l has only formal significance since operationally e z is a "third" direction, 
exclusive to both a and b. If measurements along a and b have yielded non-null polarizations, then simultaneous 
measurement along e z could only yield a null result. 

Suppose now we consider four possible polarization directions, a, a', b, and b'. Then, as is well known, the 
corresponding CHSH string of expectation values 13j, namely 



£(a, b) + £(&, b') + £(&', b) - £(a', b') , 
can be rewritten in terms of the products of our local functions A(a, /x) and B(h, /x) as 



lim 



n 

-J2 {MS) Sb(/x l ) + AJji^Bvtf) + A a ,(//)5 b (//) - A a ,{^)B h ,{^)} 



i = 1 



(34) 



(35) 



But since the standard scores ^4 a (/x) and i?b(/x) represent two independent points of S 3 , we can take them to belong 
to two disconnected "sections" of S 3 (i.e., two disconnected 2-spheres within S 3 ), satisfying 



[A n {n), B a ,(fi)] =0 V n and n' £ E 3 



(36) 



which is operationally equivalent to assuming a null result along the direction n x n' exclusive to both n and n'. If 
we now square the integrand of Eq. (|35p, use the above commutation relations, and use the fact that, by definition, 
all local functions square to unity (the algebra goes through even when the squares of the local functions are allowed 
to be —1), then the absolute value of the CHSH string leads to the following variance inequality: 



\£ (a, b) + £ (a, b') + £{a', b) - £(a', b')| 



\ 



lim 

n»l 



n 

-J2 {4+ [J^{jt% A^Qi*)][B v {fti), B h Qt*)]} 



i = 1 



(37) 



Here the classical commutators [A a (fi l ), A a >(/x 1 )] and [Bb^fi 1 ), .Bt^/x 1 )] are simply geometric measures of the 
parallelizing torsion within the 3-sphere [8j. Next, using the definitions (|25p and (PB]) for the variables A a (fi) and 
£?b (/x) and making a repeated use of the identity ( + /x-a)( + /x-a') = a • a ' — /x-(axa'), the above inequality 
can be simplified to 

\£ (a, b) + £ (a, b') + £{a', b) - £(a', b')| 



\ 



lim 



- ^ |4 + 4[ - (axa')'(b'xb) - /x 4 • 6 ] j 



<2i/l-(5xa')-(b'xb), 



(38) 



where /x l • is a null bivector. It is worth noting here that — because of the nullity of /x l • — the summation over the 
microstates /x l is not needed to arrive at the last inequality— the summation simply becomes redundant, bringing out 
the purely geometrical character of the EPR correlations [8( in addition to conforming the consistency of assumption 
(|3l>)) . Moreover, using the definitions (|I8I) for the vectors a and b, the last inequality can be further simplified to 



\£(a, (3) +£(a, (3') +£(<*', (3) ~ £(a / ,f3')\ 



^ 2 sin2(a - a') sin2(^ - J 7 ). 



(39) 
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It is now easy to see that a maximum violation of 2-\/2 of the Bell-CHSH inequality can be achieved for the angles 
(a, a', /3, /?') = (0°, 45°, 22.5°, 67.5°), precisely as observed in the experiments (HQ. On the other hand, by using 

- 1 < sin2(a - a') sin2(/3 - j3') ^ +1 , (40) 

the above inequality can be reduced to the form 

\S(a, P) + £(a, /?') + £(a', (3) - S(a', /?') | ^ 2^2, (41) 

which is precisely what is predicted by quantum mechanics, and observed in the experiments 

We are now in a position to explain where exactly the illusion of quantum nonlocality stems from The key to 
our explanation is the topology of the parallelized 3-sphere, or equivalently the algebra of rotations in the physical 
space As we saw in equation (|32j) . topologically 3-sphere is simply a 2-sphere worth of circles, but with a "twist" 



in the bundle [14| [15[ [16j ■ That is to say, as an S 1 fiber bundle over S 2 , S 3 ^ S 2 x S 1 . And this twist — analogous to 

the one in a mobius strip — is responsible for producing the right combination of polarizations, namely H — h, , H — , 

and — +, observed in the experiments To bring this out explicitly, we have modeled the physical space, not as 

H 3 , but S 3 , in terms of coordinate- free geometry of points and planes introduced by Grassmann 0]. This, in turn, 
has allowed us to introduce the raw scores observed by Alice and Bob as g/(a, fi) = (— 7-a)( + /x-a) = ±1 and 
3§(f3, fi) = (+ 1 ■ h) ( + fi ■ h) = ±1, respectively, with the random factors ( + fj, ■ a) and ( + fj, -h) representing the 
photon polarizations, and the fixed factors (— I ■ a ) and (+ 1 ■ b ) representing the analyzers used by Alice and Bob to 
measure these polarizations. Furthermore, we have taken the randomness /x — + 1 or — / shared by Alice and Bob to 
be the initial orientation (or handedness) of the entire physical space, or equivalently that of a 3-sphere. Consequently, 
once \i is given as an initial state, the polarizations along all directions chosen by Alice and Bob would have the same 
value, because /x completely fixes the sense of bivectors fi ■ n belonging to S 2 C S 3 , regardless of direction. However, 
and this is an important point, the polarization ( + fj, ■ a) observed by Alice is measured with respect to the analyzer 
(— I ■ a), whereas the polarization ( + fi ■ b) observed by Bob is measured with respect to the analyzer (+ J • b ). 

To appreciate the significance of the last point, let us begin with the scenario in which Alice and Bob are oblivious of 
each other's existence. Then, as is evident from equations (|T6]) and (fTTj) . their measurement results would depend only 
on the sense of fx. Consequently, Alice, for example, would be justified in concluding from jz/(a, fi) = + 1 that the 
local handedness associated with the points of S 3 is of counterclockwise variety, since (— 7-a)( + jU-a) = +1 holds 
for any direction a chosen by Alice if fj, = + /; and, fi, as we noted, defines a consistent sense of handedness over the 
whole of S 3 for each run of the experiment. Similarly, Bob too would be justified in concluding from /x) = — 1 

that the points of the 3-sphere are oriented in the counterclockwise sense, since (+J-b)( + /x-b) = — 1 holds for 
any direction b chosen by Bob if fj, — + 1. Furthermore, these conclusions are unproblematic in the important case 
of "perfect anti-correlation" — i.e., when b happens to be equal to a. But since we have — + 1 and 33 = — 1 in 
this case, this is, in fact, the case "H — " observed in the experiments. Analogously, for the initial orientation \x = — I 
the results of Alice and Bob would be s$ = — 1 and 38 = + 1, and that would be the case " — h" observed in the 
experiments (with equal probability). Note, however, that so far we have no indication of "nonlocality." 

But now suppose we let b — > — a. Then (+ 1 ■ b ) becomes {+ 1 ■ (— a)} = (— I ■ a), and — unlike in the previous 
case — the analyzers used by Alice and Bob become mathematically identical, with the corresponding raw scores srf 
and S% generated with the same scale of dispersion for each direction. As a result, now it would be impossible for 
Alice and Bob to observe opposite polarizations without violating the consistency of handedness defined by /x over 
the whole of S 3 . Thus, if Alice's result turns out to be (— I ■ a ) ( + ■ a ) as before, then Bob's result must also be 
(— / •a)( + /i-a), albeit for the detection of polarization along the angle a + 5 . In other words, the consistency of 
handedness throughout the 3-sphere necessitates that the only possible results Alice and Bob could observe in this 
case are either stf — + 1 and 3S — +1, or jz/ = —1 and Bd = — 1, depending on whether the initial orientation for the 

run has been /j, = + 1 or — I. These are, then, the cases "++" and " " observed in the experiments (again, with 

equal probabilities). 

Following Bell Q , such instantaneous changes in the relative measurement outcomes observed at remote stations are 



usually taken as evidence of "quantum nonlocality" 17j. In our spherical model of the physical space, however, such 
changes are determined entirely by the intrinsic geometrical and topological structures of the 3-sphere (or equivalently 
by the algebra of rotations in the physical space) without requiring any form of communication between Alice and 
Bob. To appreciate this in detail, let us express the results of Alice and Bob as limiting cases of two quaternions 
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constituting the 3-sphere: 



and 



S 3 3^(a, At) = Jim_[(-/-a)( + M-a')] = +A 



S 3 3 38(fi, fx) = Jim [(+/-b)( + /i-b')l 

b'-s-b 



(42) 



(43) 



where A = ± 1, with fi = XI, and we assume that the angles 9^ » and , between a and a ' and b and b ', respectively, 
are infinitesimally small. Moreover, as we saw in Eq. (1191) . the quaternion (+7-a)( + /Lt-a') can be decomposed as 

(+J-a)( + /i-a / ) = -A{5-5' + !• (SxS')} 

= -A{cos<9 53 , + (I-c) sin(9 55 ,}, (44) 

with c = a x a'/|a x a'|. Now suppose the pair bb' of vectors starts out being aligned with the pair aa', and then 
gets rotated counterclockwise by angle 9-^ relative to the pair a a', about the c-axis. Such a rotation within S 3 can 
be described intrinsically, as 



(+J.b)( + M-b') = {( + J-a)( + /i-a')} 



ab 

= -A{cos(0 55 



+ (J.c)sin(0- 



ab 



0>: 



(45) 



where 7?. 5 g = exp { ( Z • c ) 9~^ } = a b is the rotor that parallel transports the quaternion (+I-b)( + /x-b') from 
"point a" to "point b" within the 3-sphere QHl 



14]. But using the identity 



cos (0 s£ 



and noting that I ■ c in the limit a ' 
@(fi, /x) = litn_ [ | 

b'^b 

— > Jim_[72 sE {(+/.a)( + M-a')}] = -A 



a reduces to a null bivector, it is easy to deduce from this rotation that 
7-b)( + At-b')l = -A 



COS I 



~ £ = - A cos2(a - fi) 



(46) 



This is a probabilistic prediction that says that the relative frequency of occurring the result 38 (fi, At) = — A will vary 
with the cosine of the angle 9~^. In particular, as b = a —> — a, the result 38 (fi, A 4 ) wm change from — A to + A, and 
consequently the value of the product £/(a, fi) 38 (fi, fi) will change from —1 to +1, just as we discussed above. This, 

15[[16|. By contrast, 
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then, is a manifestation of the "twist" in the Hopf fibration of the 3-sphere mentioned above 

had we modeled the physical space as H 3 instead of S 3 , the same result would have appeared as a nonlocal effect. 

In the previous paragraphs we investigated the cases b — > a and b > — a in some detail. For the general case 
b ^ a it is convenient to view the analyzers (— I ■ 5 ) and (+ 1 ■ b ) as two different scales of measurements, as we have 
done above. Accordingly, we may think of polarization ( + fi ■ a) observed by Alice as scaled by (— I ■ a) to obtain 
the raw score srf(a, fi) = ± 1, and likewise polarization ( + fi ■ b ) observed by Bob as scaled by (+ 1 ■ h ) to obtain the 
raw score 38 (fi, fi) = ± 1. But since (— I ■ 5 ) and (+ 1 ■ h ) are two different scales measuring the same orientation of 
S 3 , the correct correlation between the raw scores can only be determined by comparing the corresponding standard 
scores, as we have done in Eqs. ([30]) to (|4"T|) . Alternatively (but equivalently) , the correct correlation between the 
raw scores can be determined by calculating their covariance divided by the product of their standard deviations, as 
explained in Eq. ([2T]) . with the product of the standard deviations in our case being (— I ■ & )(+ 1 ■ h) = ab (see also 
[18|). This product however is precisely the rotor that quantifies the twist in the fibration of S 3 , as in Eq. (|46|) above 
(see also Ref. 14{). Its value varies from +1 for b = a to —1 for b = — a and back, producing the correct combination 

of probabilities for detection of polarizations H — h, , H — , and — h observed in the experiments. Consequently, in 

accordance with our results (l33l) and (|46|) . when the raw scores srf = ± 1 and 38 — ± 1 are compared in practice by 
coincidence counts 0] 0] , the normalized expectation value of their product will inevitably yield 



£(a, fi) 



C++ (a, fi) - 


h C_ 


_(a, f3) - 


- c+- 


-(«, fi) 


- C- 


+ («, d) 


C++ (a, fi) - 


h C_ 


_(a, fi) - 


- C+- 


(a, fi) 


+ C_ 


+ («. P) 



lim 

n> 1 



n 

- V si(a, ft) 38(j3, ft 1 ) 

71 ' 



cos 2 (a 



(47) 
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where C_| (a, f3) etc. represent the number of joint occurrences of detections + 1 along a and — 1 along b etc. 

Now the important feature of the experiments performed in Orsay [4| and Innsbruck [5| was that the measurement 
settings at the remote polarizers were changed during the flight of the two particles, thus removing any possibility 
of communication between the two ends. The Innsbruck team precluded the possibility of communication between 
the polarizers by using random and ultrafast switching of the orientations of the polarizers. On each side of the 
experiment a local computer registered the polarizer orientation and the result of each measurement, with timing 
monitored by an atomic clock, and the data was gathered and compared for correlation only after the end of each run. 
These achievements prompted the principal investigator of the Orsay experiment to make the following comment: 

I suggest we take the point of view of an external observer, who collects the data from the two distant 
stations at the end of the experiment, and compares the two series of results. This is what the Innsbruck 
team has done. Looking at the data a posteriori, they found that the correlation immediately ch ang ed as 
soon as one of the polarizers was switched, without any delay allowing for signal propagation... [17| . 

But given the explanation above this immediate change in correlation should be no more puzzling than the sudden 
change in the color of Dr. Bertlmann's sock discussed by Bell In fact, in our opinion classical men like Hamilton, 
Grassmann, Clifford, and Hopf would not have been puzzled by this change at all. They would have explained the 
results observed in the experiments in exactly the same local terms as we have explained them here; namely, in terms 
of the algebra of rotations in the physical space S 3 , otherwise known as Clifford algebra Q 00101 • 
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APPENDIX: ERROR PROPAGATION WITHIN A PARALLELIZED 3-SPHERE 



In this appendix we spell out the statistical basis of the results and $Zjo§ in more detail. To this end (recalling that a 
unit 3-sphere is homeomorphic to a set of unit quaternions), consider a probability density function -P(q) : S 3 — > [0, 1] 
of random quaternions over S 3 : 

V^Mq)!! 2 I 2|Kq)|| 2 / 

It is a matter of indifference whether the distribution of q so chosen happens to be normal or not [ll|. Here 
q = s$ + w = scalar + bivector is a generic quaternion within S 3 whose mean value m(q) is defined as 

n 

m(q) = - V q< (49) 
n 

i — 1 

and whose standard deviation cr(q) (with ||q|| 2 = q^q) is defined as 



-Y, II - ™(q)ll 2 ■ (so) 

i = 1 

Needless to say, these definitions remain valid for the limiting cases of the random variable q € S 3 when 
q = € [— 1, +1] and q = w e S 2 . Now let w = p fi ■ n G S 3 be a random bivector with ||/x • n|| 2 = 1 and p £ [0, 1], 
where fi = XI is the indefinite volume form as before and n € H 3 is a unit vector. Then, using equation ([5U)) . it is 
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easy to verify that, in general, the mean m(w) of w and the standard deviation er(w) of w, respectively, would be a 
bivector and a scalar: 

m(w) = a bivector 
and ct(w) = a scalar. (51) 

Suppose now we consider another non-random bivector, v=I • n, and define a scalar number — 1 < s/ = v w < +1 , 
so that m(sf) > 0. Since v is a non-random bivector, errors generated within s/ by the random process pX stem 
entirely from the random bivector w, and propagate linearly. In other words, the standard deviations within the 
random number s/ due to the random process pX is given by 

a(sf) = vcr(w). (52) 

But since er(w) is a scalar, the typical error a(s/) generated within s/ due to the random process pX is a bivector 
[l8| . The standardized variable (which must be used to compare the raw scores s/ with other raw scores 23 ) is thus 
also a bivector: A :— s/ /cr(s/) = scalar x w. 

As straightforward as it is, the above conclusion may seem rather unusual. It is important to recall, however, that 
in geometric algebra scalars and bivectors are treated on equal footing. They both behave as real- valued c- numbers, 
albeit of different grades To appreciate the consistency and naturalness of the above conclusion, let 

s/ = /(w) = vw (53) 

be a continuous random scalar generated by the geometric product of the two bivectors v and w, as before. The 
natural question then is: How does a typical error in w governed by the probability density (|48[) — which can be 
represented by the 68% probability interval 



[to(w) — er(w), m(w) + <r(w) ] 



(54) 
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as shown in the Fig. (QJ — propagate from the random bivector w to the random scalar si , through the function 
/(w) = vw? To answer this question we note that the two end points of the interval (1541 represent two non- 
equatorial points, q~ and q + , of the parallelized 3-sphere, which is a Riemannian manifold. The geometro-algebraic 
distance between the points q and q + can therefore be defined as 

d(q~, q + ) = (q~ - q + ) x sign(q" - q+) . (55) 

Moreover, from the definition ([53)) of and a first-order Taylor expansion of the function /(w) about the point 
w = m(w) we obtain 



si = /(m(w)) + |£ 



(w - m(w)) + . . . (56) 

w — m(w) 



Now it is evident that the slope df/dw — v of this line is a constant. Therefore the mean m(&/) and the standard 
deviation of the distribution of si can be obtained by setting w = m(w) and w = er(w): 

m(rf') = /(m(w))=vm(w) = a scalar . . 

and < J (^) = faj °"( w ) =va(w) = a bivector. 

The probability distribution of srf is thus represented by the interval 

\m{st) - a(sf), m{sf) + <r{sf) } . (58) 

If we now assume that w is a unit bivector with a vanishing mean, then we have m(srf) — and <j{sd) = v, as in 
equation (|2"5)l above. 

It is instructive to note here that, geometrically, the propagation of error within S 3 is equivalent to a simple change 
in perspective: 

bivector scalar scalar bivector 

S 3 3 m(w) ± cr(w) /<W) > m(sf) ± a{sf) . (59) 

s ^ v. ^ 

quaternion quaternion 



In particular, the probability distribution of srf over S corresponding to (158)) is equivalent to that of w over S 
corresponding to (|54|) . 
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